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Abstract 

The behaviour of the mean Euler-Poincare characteristic and mean Betti's num- 
bers in the Ising model with arbitrary spin on H? as functions of the temperature 
is investigated through intensive Monte Carlo simulations. We also consider these 
quantities for each color a in the state space Sq = {—Q, —Q + 2, . . . , Q} of the 
model. We find that these topological invariants show a sharp transition at the 
critical point. 

Keywords: Phase transitions and critical phenomena. Topological invariants (Euler- 
Poincare characteristic, Betti's numbers). 

1 Introduction 

Knowledge about the spatial structure of systems is a subject that attracts more and 
more interest in statistical physics. For example, the study of morphological features has 
proved very useful to identify and distinguish various formation processes in contexts as 
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varied as the study of large scale distributions of galaxies, or the investigation of micro- 
emulsions [TJ. Other examples of situations where the focus on the spatial structure is of 
primary importance to understand the physical properties of a given system are provided 
by dissipative structures in hydrodynamics, Turing patterns occurring in chemical reac- 
tions, or transport properties of fluids 2\. For instance, the transport of a fluid over a 
porous substrate depends crucially on the existence of a connected cluster of pores that 
spans through the whole system [3j. Gaining insight into the spatial structure of the most 
commonly used lattice-models in statistical physics -i.e. the Potts, Clock and generalized 
Ising models- is also a matter of interest in the field of image-analysis, which makes ex- 
tensive use of these models to simulate noise and clean dirty images jlj. In this context, 
mathematical quantities that characterize the typical spatial structure of the model seem 
to be the most natural tuning parameters as opposed to the usual thermodynamic quan- 
tities. 

Consequently, efforts have recently been made [3 El El to study the typical behavior of 
Minkowsky-functionals (or curvature integrals) which quantify the geometrical properties 
of a given system. Particular interest has been devoted to the study of the Euler-Poincare 
characteristic which turns out to be a relevant quantity in various models : let us indeed 
recall, that for the problem of bond percolation on regular lattices, Sykes and Essam 
[7] were able to show, using standard planarity arguments, that for the case of self-dual 
lattices (e.g. Z 2 ), the mean value of the Euler-Poincare characteristic changes sign at the 
critical point, see also jH]. 

More recently Blanchard, et al. showed jU] using duality and perturbative arguments, 
that for the two dimensional random cluster model (the Fortuin-Kasteleyn representation 
of the q-state Potts model) the mean local Euler-Poincare characteristic is either zero or 
exhibits a jump at the self-dual temperature. More precisely, it vanishes when q — 1 and 
exhibits a jump of order 1/y/q when q is large enough. Similar results were shown in 
higher dimensions d. 

This paper aims to investigate the spatial structure of the Ising model with arbitrary 
spins on Z 2 which, unlike the Potts-model, does not have a color-symmetry. We intro- 
duce colored complexes (i.e. the sets of sites for a given configuration that belong to the 
same color) and study the behavior of the associated local mean Euler-Poincare charac- 
teristics. As a result of our numerical computations, we obtain that the mean value of 
the Euler characteristic per site vanishes below the critical temperature and is positive 
above this temperature. As we shall see, this behaviour can be related to the behaviours 
of topological invariants, the Betti's numbers. 

The paper is organized as follows. In Section 2, we introduce the model and give the defi- 
nitions of the principal quantities of interest. Section 3 is devoted to the results obtained 
by intensive Monte-Carlo simulations. Concluding remarks are given in section 4. 



2 Definitions 

To introduce the generalized Ising model considered by Griffiths' in :9j, we associate to 
each lattice site x G Z 2 a spin a x in the set Sq = {—Q, —Q + 2, . . . , Q}, Q = 1,2,..., 
(with cardinality \Sq\ = Q + 1). The Hamiltonian in a finite box A C Z 2 is given by: 




(1) 
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where J is a positive constant that we will take equal to 1 in the next section, and the 
sum runs over nearest neighbors pairs of A. The associated partition function is defined 
by: 

Z A (P)= J2 e ~" HA ( 2 ) 
Ox- € Sq 

x e A 

where the sum is over all configurations cr A = {a x } xEA and we let 

e Sq 

x e A 

denote the (corresponding expectation) of a measurable function /. Among the quantity 
of interest, we shall first consider the magnetization 

Mfl = |^<X>»>a(/J) (4) 

where |A| is the number of sites of A. Let us recall that this system exhibits a spontaneous 
magnetization at all inverse temperatures greater than the critical temperature of the Ising 
model (Q — 1), as proved in 0. 

We next introduce others quantities, namely the mean Euler characteristic per site. 
Consider a configuration <x A , and let L(<x A ) be the set of bonds (unit segments) with 
endpoints x, y, such that a x = a y , and -P(cr A ) be the set of plaquettes p with corner 
x,y,z,t, such that a x = o y = o z = o t . To a given configuration cr A , we associate (in 
a unique way) the two-dimensional cell-sub complex C(cr A ) = {A, L(cr A ), P(cr A )} of the 
complex C = {A, L A , P A } where L A is the set of bonds with both endpoints in A and P A 
is the set of plaquettes with corners in A. The Euler characteristic of this subcomplex is 
defined by: 

X {<T A ) = \A\-\L( ( r A )\ + \P{ t T A )\ (5) 

where \E\ denotes hereafter the cardinality of the set E. It satisfies the Euler-Poincare 
formula 

X (<T A ) = 7r°(0- A ) - Tz\cT h ) + 7T 2 (0- A ) (6) 

where 7r°(cr A ) and 7r 1 (cr A ) are respectively the number of connected components and the 
number of independent 1-cycles of the subcomplex C(cr A ). Here 7r 2 (cr A ) = because 
C(cr A ) has no 2-cycles. 
Notice that 

xOa) = X^>i°c(°"a) ( 7 ) 

xgA 

where the local Euler characteristic X Xi i oc (cr A ) is given by: 

1 . 1 - . , 

X*,ioc(°a) = i - - Yl 5 ^ x ' °y} + 4 II az ) ( 8 ) 

Here, V(x) is the set of nearest neighbours of x, the second sums is over plaquettes 
containing i as a corner, and 5 is the kronecker symbol: 5(o~ x ,a y ) — 1 if a x — o~ y , and 
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^{^xi^y) = otherwise. We define the mean Euler characteristic per site and the mean 
Betti's numbers per site by: 

XtSP) = i|t(x^a))a(/3)(1) 
= |^<7rV A )) A (/3)(2) 
= ^(^(^a))a(/?)(3) 
The following limits exist and coincide: 

lim xM = lim (x 0) ioc(o-a))a(/3) (9) 

Here the first limit is taken in such a way that the number of sites of the boundary of A 
divided by the number of sites of A tends to 0, while there is no restriction in the second 
limit. 

The existence of these limits are simple consequences of Griffith's inequality JU] which 
states that: 

(«)a(/3) > «)a(/3)«)a(/3) (10) 
This implies the following monotonicity properties with respect to the volume: 

(<t>™)a(/3) < (<#r?>A'(0)if A c A' 
These properties give the existence of the infinite volume limit 

(«)a(/3) := lim 2 «<>a(/?) 
w a^z 2 m 

through sequences of increasing volumes. They imply also, following ^U|, that the limiting 
state is translation invariant. It is then sufficient to observe that: 

N Ubes Q ,bM b - ^ - a v^ 
o{a x ,a y ) = — • 



Ylbes Q ,b^o b 

and thus can be written as a (finite) linear combination m€S Q -Vfc "™ "™ with finite 
coefficients \ a ,b to conclude the existence of the second limit in (J0J), since also by (JHJ) 
the local Euler characteristic is a (finite) linear combination of the <J™a™ with finite 
coefficients. The translation invariance allows finally to show that this second limit in © 
coincides with the first one under the above mentioned conditions. 

We shall also consider the subcomplexes of C(cr\) attached to each color. Namely, for 
a given configuration a a and a given color a G Sq, we let A a be the subset of the box A 
such that a x = a, L a (a\) be the subset of bonds of L(a\) with endpoints x, y, such that 
o'x = cr y = a, and P a (er/C) be the subset of plaquettes of P(<ta) with corner x, y, z, t, such 
that a x = o y = a z = <j t = a. For each color a, we associate to a given configuration a\, the 
two-dimensional cell-subcomplex C a (cr A ) = {A a , L a (cr A ), P a (cr A )}. Obviously C a (cr A ) is 
a subcomplex of C(cr A ) and hence a subcomplex of C and C(cr A ) = U ae s Q C a (cr A ). 

The Euler characteristics of these subcomplexes are defined by: 

X a (<T A ) = \A a \ - \L a (<T A )\ + \P a (<r A )\ (11) 
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They satisfy the Euler-Poincare formula 

X (<r A ) = tt 'Va) - tt^^o-a) + tt 2 '> a ) (12) 

where 7r°' a (cr A ) and 7r 1,a (cr A ) are respectively the number of connected components and 
the maximal number of independent 1-cycles of the subcomplex C a (cr A ). Here again 
7r 2 ' a (cr A ) = because C a (cr A ) has no 2-cycles. 
Note that 

As for the total Euler characteristic, the Euler characteristic of each color can be written 
as a sum over local variables: 

X>a) = $>^°c( ct a) ( 14 ) 

x€A 

where 

X^IocOa) = S(a x ,a) - - ^ 6(a x ,a)5(a y ,a) + ^^2Y[6(<r y ,a)5(<T g ,a) (15) 

y£V{x) p3x yz&p 

We define the mean Euler characteristic per site and the mean Betti's numbers per 
site of the colored subcomplexes: 

XX(P) = ^<X>a)>a(/3)(4) 

^{P) = ^°'Va))a(/3)(5) 

vrI' a (/3) = ^(^Va)^)^) 

One can show that the following limit exist and coincide: 

Jim X %(J3) = lim 2 (xo i1oc (ct a )) a ((3) (16) 

by using the relation 

, - Hb£S Q ,b^a( a x ~ b) 

5{<7x,a =-fr 1 K\ 17 

and arguing as before. 



3 Results 

We have performed intensive Monte Carlo simulations to study the thermodynamic be- 
haviour of this model. As expected, we find that it belongs to the 2D-Ising universality 
class for all tested values of Q with the same critical exponents. The inverse critical tem- 
perature (3 C at which second order phase transition occurs for Q ranging from 2 upward 
tends to with increasing values of Q. 

Using Binder cumulants [II] (see also P2J, d )> we find (3f =1 = 0.4405 ± 0.0004, 
pQ= 2 = 0.1475 ± 0.0005, p®= 5 = 0.0316 ± 0.0005. The critical exponent for the magneti- 
zation is found equal to 0.125 ± 0.001 and the one for the susceptibility to 1.76 ± 0.002. 
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Figs. 1 and 2 show the finite size scaling plots of the magnetization and the suscepti- 
bility for Q = 2 and Q = 5 for different system sizes L: |A| = L x L. 
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Figure 1: Scaling plots of the magnetization for Q = 2 and Q = 5. (M) L 1 ^ 8 as function 
of (3*L. (3* = is the reduced inverse temperature. 
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Figure 2: Scaling plots of the susceptibility [M — (M)] 2 x L 7 / 4 as function of for 
Q = 2 and Q = 5 for different system sizes. 



In the following we present computed mean values of the Euler characteristic and 
related Betti's numbers per site and per color for several values of Q as function of the 
inverse temperature j3. 

Periodic boundary conditions have been used which obviously lead to a symmetry 
between colors a and— a. The equilibrium configuration presents a preference for positively 
or negatively oriented spins. We have chosen in the statistics, the first one, i.e. the one 
corresponding to +Q boundary conditions. 

Figure El shows the behaviour of the mean values of the global Euler characteristic Xa 
and Betti's numbers 7r° and tt\ as function of j3 for Q = 1,2,5 and |A| = 100 2 . Recall 
from (JHJ) that 7r° (/?) is nothing but the number of connected components of the complex 
and tt\(/3) is the number of independent cycles, i.e. the number of 1— dimensional holes 
in the complex. 
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Seemingly, XaW) appears as a decreasing function of (3 and vanishes for (3 > (3 C . The 
decrease of Xa(P) f° r /3 < /?c is associated both to a decrease of 7T A (/3) and an increase of 
7r\(P). In addition one observes that for (3 > (3 C these two quantities decrease and seem 
to compensate each other. This would explain the vanishing of Xa(P) f° r P ^ Pc- As 
we shall see this compensation can be understood intuitively by looking at the behaviour 
of the colored Euler-Poincare characteristic Xa(P) an d colored Betti's numbers 7r A ' a (/?), 
vr A ' a (/3), a G S Q , Q = 1,2,5. 



<:-./ 




0,1 0,2 0,3 0,4 0,5 0,6 0.7 0.3 Q Q2 0,4 0,6 0.8 1 

b b 



Figure 3: Mean values of the global Euler characteristic x an d global Betti's numbers Ilo 
and 111 for Q — 1, 2, 5 and lattice size |A| = 100 2 . 



Let us look at the Ising case Q — 1 for which these quantities are presented in Figure EJ 
We observe that for (3 < f3 c the quantities x\{0)i 7r A <I (/^)' ""/(^l do no ^ depend on colors 
a = ±1. However they exhibit a sharp transition in the color variable right at the critical 
point. 

It is also seen that above (3 > (3 C , 7r A +1 ~ 7r A ' _ ~ and 7r A ' +1 ~ tt a ~ . Therefore 7r A ~ 
ti\ which implies that in this range of temperature the global Euler characteristic vanishes. 
Alternatively the behaviour of the colored Betti's numbers explain why XaH^) — ~\ r A^P) 
which, in view of formula f|13j) . gives also that the global Euler characteristic vanishes. 
These facts can be understood intuitively at low temperatures if one has in mind the usual 
picture of islands of (— ) inside a sea of (+). Indeed, for such configurations, say a a, one 
has 7r A ' +1 (cr A ) ~ 7T A '~ (<x A ) — 1 and each island gives a contribution 1 to both 7r A +1 (cr A ) 
and 7r A _1 (cr A ). 
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On the other hand, when (3 < (3 C there exists only one (disordered) phase with inde- 
pendent colors. This corroborates the fact that X%{0)i Ta"^) anc ^ ^iCiP) do not rely 
on colors in this domain of temperature. In addition, the disorder increases when beta 
decreases, and in this domain of temperature it is natural to expect that the number of 
connected components of each color increases, while the number of independent cycles (of 
each color) decreases. 
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Figure 4: Mean values of the colored Euler characteristic x a {/3) and Betti's numbers 
Tc°' a ((3) and tt 1 '"^) for Q = 1 and |A| = 100 2 . 



Figure El shows the results for the case Q = 2 (a = 0, ±2). Considering first the colors 
a = ±2, one can formulate the same reasoning as for the case Q = 1. The quantities of 
interest behave similarly as before except for the amplitudes. 

The color a = does not present any particular feature. 7r A ' (/?) ~ translates the 



fact that this color is dominated by the two extremal states a 
13. This is reflected by noting that Xa(0) — 



-2, +2 for all values of 



~ tt\'"(/3) in Figure El 
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Figure 5: Mean values of the colored Euler characteristic x a {(3) and Betti's numbers 
7T°' a ((3) and 7T 1,a G0) for Q = 2 and |A| = 100 2 . 



In Figure |H1 we present our results for the case Q = 5 (a = ±5, ±3, ±1). We observe 
that for p > (3 C the main contribution to the number of connected components comes 
from color a = —3. 
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Figure 6: Mean values of the colored Euler characteristic x a {(3) and Betti's numbers 
TT°' a ((3) and tt 1 '"^) for Q = 5 and |A| = 100 2 . 



Figure U\ below gives an enlargement of Figure El showing the behavior of the other 
colors. Here also a sharp transition is observed at the critical point. 
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Figure 7: Enlargement of Figure El near the critical point. 

All numerics have been performed with careful thermalization control (up to 5.10 6 
MC steps before measurements). Data simulation errors were controlled through bin- 
ning analysis. In all figures, the error bars are smaller than the sizes of the symbols. 
The Hoshen-Kopelman cluster search algorithm ^Hj for colored connected components 
identification has been used. 

4 Conclusion 

In |H] we have already shown that for the Random Cluster Model, the mean value of the 
Euler-Poincare characteristic with respect to the Fortuin-Kasteleyn measure exhibits a 
non trivial behaviour at the critical point. Namely, there we found that, in dimension 2, 
it is either zero or discontinuous at the transition point. 

Here, for the Ising model with arbitrary spins on Z 2 , our numerics show that the 
topology of configurations reveals the signature of the phase transition. In particular this 
leads us to following conjecture 

Conjecture: in the thermodynamic limit, Xatz 2 (<^) > for /3 < /5 c and X\^ 2 ^) = ® 
lor 1> i. 

We think that the vanishing of X^^iP) can be shown at large (3 by perturbative 
Peierls type arguments in view of the expression (|T5|) of the local Euler characteristic. 

It would be interesting to study also these quantities in models with color symmetry 
like the Potts or Clock models [TH] . 
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